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ABSTRACT 

We examine the selection characteristics of infrared and sub-mm surveys with 
IRAS, Spitzer, BLAST, Herschel and SCUBA and identify the range of dust tem- 
peratures these surveys are sensitive to, for galaxies in the ULIRG luminosity range 
(12< log (Ltr/Lq) <13), between z=0 and z=4. We find that the extent of the red- 
shift range over which surveys are unbiased is a function of the wavelength of selection, 
flux density limit and ULIRG luminosity. Short wavelength (A <200 /im) surveys with 
IRAS, Spitzer /MIPS and Herschel /PACS are sensitive to all SED types in a large 
temperature interval (17-87 K), over a substantial fraction of their accessible redshift 
range. On the other hand, long wavelength (A>200/im) surveys with BLAST, Her- 
schel /SPIRE and SCUBA are significantly more sensitive to cold ULIRGs, disfavour- 
ing warmer SEDs even at low redshifts. In order to evaluate observations in the context 
of survey selection effects, we examine the temperature distribution of 4 ULIRG sam- 
ples from IRAS, Spitzer /MIPS, BLAST and SCUBA. We find that the lack of cold 
ULIRGs in the local (z^O.l) Universe is not a consequence of the selection and that 
the range of ULIRG temperatures seen locally is only a subset of a much larger range 
which exists at high redshift. We demonstrate that the local luminosity-temperature 
(L-T) relation, which indicates that more luminous sources are also hotter, is not 
applicable in the distant Universe when extrapolated to the ULIRG regime, because 
the scatter in observed temperatures is too large. Finally, we show that the difference 
between the ULIRG temperature distributions locally and at high redshift is not the 
result of galaxies becoming colder due to an L-T relation which evolves as a function of 
redshift. Instead, they are consistent with a picture where the evolution of the infrared 
luminosity function is temperature dependent, i.e. cold galaxies evolve at a faster rate 
than their warm counterparts. 

Key words: galaxies: general galaxies: high-redshift galaxies: starburst galaxies: 
evolution infrared: galaxies 



1 INTRODUCTION 

Infrared (IR) radiation was first associated with individual 
galaxies in the late 1960s (e.g. Johnson 1966; Low & Tucker 
1968; Kleinmann & Low 1970), laying out the path for sub- 
sequent observations by the first infrared space observatory, 
the InfraRed Astronomical Satellite (IRAS). The IRAS all 
sky survey (Soifer, Neugebauer & Houck 1987) established 
the existence of a class of galaxies whose bolometric energy 
output emerges almost entirely in the infrared (Soifer et al. 
1984a; Sanders & Mirabel 1996) and amongst those, a pop- 
ulation of ultraluminous infrared galaxies (ULIRGs; Soifer 
et al. 1986; 1987) with total infrared luminosities (Lir, 8- 
1000 fim) greater than 10 12 L Q . 

ULIRGs are amongst the most actively star-forming 
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galaxies in the Universe, with a non-trivial fraction of them 
also hosting AGN (>30 per cent, e.g. Franceschini et al. 
2003, Kartaltepe et al. 2010, Symeonidis et al. 2010). Stud- 
ies have shown that although ULIRGs are rare in the local 
Universe (e.g. Kim & Sanders 1998), they become more com- 
mon at high redshift, concomitant with the evolution of the 
infrared luminosity function and the increase in the total in- 
frared luminosity density with redshift (e.g. Takeuchi, Buat 
& Burgarella 2005). Sub-mm and radio surveys show that 
ULIRGs are already in place at z~2-4 (e.g. Kovacs et al. 
2006; Coppin et al. 2008; Seymour et al. 2008) and that to- 
gether with the luminous infrared galaxy (LIRG, Lir=10 11 - 
10 12 Lq) population, they dominate the cosmic infrared lu- 
minosity density at z>l (e.g. Lagache et al. 2004; Perez- 
Gonzalez et al. 2005). 

The properties of dust in ULIRGs has been the sub- 
ject of numerous investigations because, apart from being 
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responsible for most of the energetic output, dust is both 
a by-product and catalyst for star-formation. The temper- 
ature dust can reach in different parts of the interstellar 
medium (ISM) is related to its distribution, abundance and 
type, as well as the star-formation rate and efficiency. Initial 
studies showed that ULIRGs were characterised by warm av- 
erage dust temperatures (e.g. Soifer et al. 1984; Klaas et al. 
1997). However, these results were based on local ULIRGs 
only and early observations with the Submillimeter Com- 
mon User Bolometer Array (SCUBA; Holland et al. 1999) 
revealed a different picture of the distant Universe (e.g. 
Hughes et al. 1998; Eales et al. 1999, 2000). One of the most 
important outcomes of SCUBA blank-field surveys was dis- 
covering the existence of ULIRGs characterised by colder 
average dust temperatures than those seen locally (Kovacs 
et al. 2006; Pope et al. 2006; Huynh et al. 2007a; Coppin et 
al. 2008). 

The abundance of cold ULIRGs at high redshifts has 
been commonly ascribed to the predisposition of SCUBA 
surveys to detect such objects. Analogously, our view of lo- 
cal ULIRGs has been influenced by two factors, one of which 
relates to the potential selection effects introduced by IRAS 
surveys and the other to the apparent correlation between 
infrared luminosity and dust temperature. Analysis of lo- 
cal spectral energy distributions (SEDs) has revealed that 
more luminous galaxies are warmer than their lower lumi- 
nosity counterparts, pointing to a luminosity-temperature 
(L-T) relation (e.g. Dunne et al. 2000; Dale et al. 2001; Dale 
& Helou 2002; Chapman et al. 2003), which various works 
have incorporated in their analysis of the local luminosity 
function (e.g. Chapman et al. 2003; Blain, Barnard & Chap- 
man 2003; Lewis, Chapman & Helou 2005; Chapin, Hughes 
& Aretxaga 2009). At the same time however, it has been 
suggested that the lack of cold ULIRGs in the local Universe 
could be due to a selection effect, whereby the IRAS 60/im 
waveband is expected to favour the selection of sources with 
warmer dust temperatures (e.g. Blain et al. 2004a). 

Recently some progress in resolving these issues has 
been made. Clements, Dunne & Eales (2010) demonstrated 
that the submm properties of local ULIRGs are consis- 
tent with previous results which described them as having 
warmer temperatures than sub- mm galaxies (SMGs) of com- 
parable luminosity. In addition, Symeonidis et al. (2009) and 
Seymour et al. (2010) showed that cold ULIRGs are not 
unique to the redshifts probed by SCUBA; they are already 
present at z< 1 and seem to be the main contributors in the 
strong evolution of the ULIRG energy density with redshift. 
Nevertheless, until now, the variation in the sensitivity of dif- 
ferent instruments has meant that ULIRG samples have not 
been uniformly mapped out in redshift. Herschel (Pilbratt 
et al. 2010) is expected to enable a more quantitative com- 
parison than previously possible between local and distant 
ULIRGs, but such an endeavour will require a fundamental 
understanding of a) the range of ULIRG SED types and b) 
what slice of the ULIRG population each survey can probe. 

In this paper we address these topics by firstly focus- 
ing on the IRAS view of the local Universe (Section [2]): the 
temperature range that IRAS surveys are sensitive to and 
the SED types of local ULIRGs (characterised by their peak 
wavelength, A pca k). We subsequently examine what types of 
ULIRG SEDs are traced by Spitzer, BLAST and SCUBA at 
z ^ 4 and investigate the differences between the local and 



high redshift samples (Section [3]). In Section [4] we examine 
the potential of Herschel and SCUBA 2 surveys to charac- 
terise the ULIRG population in the distant Universe and 
finally, our discussion and conclusions are presented in Sec- 
tions [5] and [6] Note that in order to discuss the differences 
in Ap Ca k distributions between the samples we examine in 
this paper, we occasionally adopt a definition of warm and 
cold SEDs as those that peak shortwards and longwards of 
90 /im respectively. Although, this division is not used as 
a substitute for a more quantifiable characterisation of the 
populations, it allows us to place our results in the context 
of earlier work, where the warm and cold terminology (or 
variations of it, such as "cool" or "hot") is used extensively 
to qualitatively describe infrared SEDs (e.g. Dunne et al. 
2000; Chapman et al. 2002, 2003; Xilouris et al. 2004; Casey 
et al. 2009). Our choice of 90 /im corresponds to a greybody 
temperature of 29 K which is (i) roughly compatible with 
the transition wavelength between IR emission from warm 
dust in young star-forming regions and a cooler component 
associated with more extended dust heated by the interstel- 
lar radiation field (Lonsdale-Persson & Helou 1987; see also 
Seymour et al. 2010); (ii) consistent with previous studies of 
local ULIRGs where the cold dust component is placed be- 
low 30 K (Klaas et al. 2001, Spinoglio, Andreani & Malkan 
2002) and 28 K (Hawkins et al. 2001) and the warm compo- 
nent above 30 K (Vlahakis, Dunne & Eale 2005); (iii) in the 
middle of the range of temperatures probed by "cold" high- 
z SCUBA SMGs and "warm" local (IRAS) ULIRGs (e.g. 
Clements, Dunne & Eales 2010) and (iv) about halfway in 
the range of SEDs of the Symeonidis et al. (2009) 70 (im 
sample. 

Throughout we employ a concordance consmology of 
H =70kms" 1 Mpc _1 , fi M =l-fU=0.3. 



2 LOCAL ULIRGS AS SELECTED BY IRAS 

In Symeonidis et al. (2009, hereafter S09) we examined the 
SEDs of local (zSCO.l) galaxies selected at 60 /xm down to 
5.24 Jy from the IRAS Bright Galaxy Sample (BGS, Sanders 
et al. 2003). Their rest frame, [vh v ] SED peak wavelength 
(Apeak) and total infrared luminosity (Lir, 8-1000 /im) were 
estimated by fitting available photometry with the Sieben- 
morgen & Kriigel (2007, hereafter SK07) SED templates, 
adopting A poa k as a proxy for dust temperature. We found 
that for a given infrared luminosity, local ULIRGs showed 
warmer SEDs which peak at shorter wavelengths than what 
was observed for a sample of higher redshift (0.1<z<1.2) 
Spttzer /MIPS 70 /im-selected sources. Here we re-examine 
selection issues and opt to provide a more stringent analysis 
of local galaxy SEDs (i) by investigating the selection char- 
acteristics of IRAS surveys and (ii) by extending our work 
to much fainter IRAS fluxes. 

We examine the sensitivity of 60 /im IRAS surveys to 
various SED types as a function of redshift, for a given flux 
density limit. SEDs are characterised by their peak wave- 
length, however to facilitate interpretation of our results, 
we relate A pe ak to temperature, using the Wien displace- 
ment law for a vi v grey body of dust emissivity (9=1.5: 

be 

T(K) ~ — (1) 

1 ' (4 + P)xkX Apeak ^ ' 

where h is the Planck constant, c is the speed of light in a 
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Figure 1. Plots of SED peak wavelength (A poa k) (and temperature, right-hand y-axis) vs redshift for surveys with IRAS for three 60 ^m 
flux density limits, 5.24, 1.2 and 0.36 Jy, corresponding to the BGS, the Fisher et al. (1995) sample and the IIFSCz sample respectively. 
Left panels and right panels correspond to logLijj.=12 and logLm,=12.5 ULIRGs respectively. Thick red patterns designate the region 
where all templates with a certain Apeak are detected at a given redshift, whereas unshaded areas correspond regions where no templates 
are detected. Blue dashed regions indicate that for a given A pea k only a fraction of the SK07 templates are recovered. 



vacuum and k is the Boltzmann constant. We choose 60/im 
flux density limits of 5.24, 1.2 and 0.36 Jy, corresponding to 
the BGS, the Fisher et al. (1995) sample and the Imperial 
IRAS-FSC Redshift Catalogue (IIFSCz; see below) respec- 
tively. For our purposes, we use the entire SK07 library, com- 
posed of 7088 SED templates which peak in the 9.4-175 fim 
range. Many templates are characterised by more than one 
dust temperature component and hence the SK07 library is 
able to match a large variety of systems. We characterise the 
SEDs via their peak wavelength which is a unique feature 
and represents the region where the bulk of the emission 
originates. Figure [1] is created as follows: (i) all templates 
are normalised to log (Lir/Lq, 8-1000 fim) = 12 and 12.5, (ii) 
they are redshifted to z=0.05 up to z=0.8 in steps of 0.05 
and (iii) the observed 60 /jm flux density is estimated for 
all templates at all redshifts. The A pca k-z parameter space 
in figure [T] illustrates the fraction of templates with a cer- 
tain Ap Ca k that are detected at a given redshift. Red thick 
patterns mark the regions where all templates are detected, 
whereas the blue dashed pattern indicates intermediate re- 
gions where only a fraction of the SK07 templates are re- 



covered. In terms of the latter, the detection rate relates 
to variations in SED shape: when 60 /im probes the Wien 
side of the SED, only SEDs with a shallow mid-IR slope 
are likely to be detected close to the surveys' flux density 
limit, whereas when it probes the Rayleigh- Jeans side, only 
templates with a substantial amount of flux after the peak 
will be recovered. In the unshaded areas, no templates are 
detected for a given flux density limit and ULIRG luminos- 
ity. For our purposes, we classify a survey as unbiased and 
complete if it can detect all SEDs between 30-150 /im (17- 
87K). 

Figure [T] demonstrates that the two deeper IRAS sur- 
veys are unbiased over certain redshift ranges. Depend- 
ing on the ULIRG luminosity, this extends to z=0.05 and 
z=0.1 for the 1.2 Jy survey and to z=0.1 and z=0.15 for 
the 0.36 Jy survey, although for the latter, ULIRGs with 
SED peak wavelengths in the 30-130 fim range can be de- 
tected up to z—0.2. For the shallower BGS, the detection of 
log (Ltr/L0)~12 ULIRGs is only complete at hotter tem- 
peratures (Apeak =40-60 /im), but the same survey is un- 
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biased up to z=0.05 for sources in the more luminous, 
log (Li R /L Q )~12.5, range. 

In the context of the selection function in Fig. [1] we 
examine the SEDs of IRAS-selected ULIRGs, using the 
Imperial IRAS-FSC Redshift Catalogue (IIFSCz, Wang & 
Rowan-Robinson 2009), which is a sub-sample of the IRAS 
Faint Source Catalogue (FSC, Moshir 1992) and about 14 
times deeper than the BGS. We keep all sources down to 
a 60 /im flux density limit of 0.36 Jy, where the sample is 
90 per cent complete with 90 per cent of the sources hav- 
ing either spectroscopic or photometric redshifts. For SED 
analysis, we also require 25 and 100 /jm detections and cal- 
culate Lir and A pea k by fitting the SK07 models as outlined 
in S09, after colour-correcting the photometry as specified in 
the IRAS Explanatory Supplement (Beichman et al. 1988). 
The effects of the additional 25 and 100 pm selection are as- 
sessed as follows: using the 60 /jm flux density we calculate 
Lir according to Lir=1.7x AL/60 (Chapman et al. 2000) and 
examine the 25 and 100 /im detection rate of sources in the 
ULIRG (12< log (Lm/L©) <13) luminosity range. We find 
that z<0.1 ULIRGs are detected in all 3 bands, however for 
z>0.1 sources the detection rate is 75 per cent at 100 /mi, 
with 25 per cent of those also having 25 /im data. As a result 
and in order to work with a complete sample, we focus all 
subsequent analysis on local (z^O.l) ULIRGs only, leaving 
us with 27 sources — cf. 15 local ULIRGs in the BGS in 
S09. 

Figure shows the rest-frame SED peak wavelength 
(grey body temperature on the right-hand y-axis) against 
total infrared luminosity (Lir, 8-1000 /mi) for the 27 local 
IIFSCz ULIRGs. For comparison we also include the local 
luminosity-colour relation from Chapin, Hughes & Aretxaga 
(2009, hereafter CHA09). Colour, defined as log [f 60 /fioo], is 
equivalent to temperature, so this relation describes sources 
as being hotter if more luminous; see Appendix |X] for details 
on converting colour to SED peak wavelength (and temper- 
ature) . 

The IIFSCz comprises about 50 per cent more ULIRGs 
at z<0.1 than the S09 BGS, including some sources with 
Apeak > 70 /mi, however, the majority (80 per cent) are con- 
sistent with the temperature distribution we saw in S09. In 
addition, 70 per cent are within la of the local luminosity- 
temperature (L-T) relation of CHA09. Despite the fact that 
the 0.36 Jy survey is sensitive to all ULIRG SED types 
at z^O.l (Fig. [TJ, the average peak wavelength is short: 
(Apeak) =62±2 fim with an intrinsic dispersion <7=12/mi; 
(T)=42±2K with a~8K. In addition, A pC ak extends only 
up to 90 /im (albeit large errors on some sources). Our re- 
sults confirm that cold ULIRGs are rare in the local Universe 
(see section Q] for our definition of cold). 



3 ULIRGS AT HIGH REDSHIFT 

3.1 The selection function of Spitzer, BLAST and 
SCUBA 

To quantitatively describe the selection issues relating to 
various infrared and submm surveys, we develop a schematic 
representation of the selection functions of Balloon-borne 
Large Aperture Submm Telescope (BLAST; Pascale et al. 
2008), SCUBA and Spitzer 70 and 160 pm surveys for 




Figure 2. SED peak wavelength (and dust temperature derived 
using equation [l] on the right-hand y-axis) versus total infrared 
luminosity (8-1000 /im) for the 27 local (z<0.1) ULIRGs selected 
at 60 /im (f 60 >0.36Jy) from the Imperial IRAS-FSC Redshift 
catalogue. The solid and dashed lines are the local luminosity- 
colour relation and lcr limits from Chapin, Hughes & Aretxaga 
(2009), converted into a luminosity- A pca k relation using equation 
lAll in Appendix |X] Note that cold ULIRGs are rare in the local 
Universe. 



log (Lir/L q ) = 12.5 ULIRGs up to z=4. Figure |3] is built 
following the method outlined in section [5] where the A pea k- 
z parameter space illustrates the temperatures these sur- 
veys are sensitive to at any given redshift and for a given 
flux density limit. As in Fig. [1] unshaded regions and those 
shaded red correspond to no and all SK07 templates be- 
ing detected respectively. In the blue dashed areas, only 
some of the templates are recovered. For the submm sur- 
veys (BLAST and SCUBA), this corresponds to the detec- 
tion of templates with significant flux in the wavelength 
range probed by these instruments, which is always long- 
ward of the peak. For Spitzer 70 and 160 pm the same ap- 
plies only if the SED peaks shortwards of the wavelength of 
selection, otherwise it is the mid-IR continuum slope that 
determines the detection rate (see also section 0. The flux 
density limits that we use for each survey are: lOmJy for 
MIPS 70 pm, matching our observations in S09, 50mJy for 
MIPS 160 /im (Dole et al. 2004), 40, 30 and 20mJy for 
the 250, 350 and 500 /im BLAST bands respectively (Dye 
et al. 2009, hereafter D09) and 2mJy for SCUBA (Cop- 
pin et al. 2008, hereafter C08). In Fig. we also overlay 
data for log (Lir/Lq) sJ12.5 ULIRGs from the S09 70 pm 
sample, SCUBA measurements from Coppin et al. (2008) 
and BLAST measurements from Dye et al. (2009). For the 
SCUBA and BLAST data we use equation[T]to convert their 
derived temperatures to A pca k. Ideally, we should be look- 
ing at ULIRGs in a narrow slice around log (Lir/Lq) =12.5 
to be consistent with the selection function displayed, how- 
ever the samples are too small and hence we plot all 
12^ log (Lir/L0)^12.5 sources. In addition, all subsequent 
estimates of the average peak wavelength and temperature 
refer to these log (Lir/Lq)^12.5 sources only, in order to 
achieve a fair comparison with local ULIRGs whose lumi- 
nosities do not extend beyond log (Lir/L0)=12.5. 
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Figure 3. Plots of SED peak wavelength (and grey body temperature on the right-hand y-axis) vs redshift for surveys with Spitzer, 
BLAST and SCUBA for log (L m /L )=12.5 ULIRGs, using flux density limits of 10, 50, 40, 30, 20 and 2mjy for the 70, 160, 250, 350, 
500 and 850 fim bands respectively. We also plot all logLiR ^12.5 ULIRGs from the samples of Symeonidis et al. (2009) (70 and 160 fim 
panels), Dye et al. (2009) (250, 350 and 500 fim panels) and Coppin et al. (2008) (850 fim panel). In the red thickly patterned regions 
all templates with a particular Ap ea k are detected at a given redshift, whereas white areas indicate that no templates are detected. Blue 
dashed regions indicate that only some of the SK07 templates with a particular A pea k are detected at that redshift. 



The MIPS 70 fim survey can potentially detect objects 
up to z=2, but with severe incompleteness above z=l — it is 
unbiased up to z ~0.7 and sensitive to warmer SEDs above 
that redshift. Our predictions are consistent with results 
from Casey et al. (2009) who find all z>l, 70 /im-detected 
ULIRGs to be hot. Surveys at 160 fim have a similar selec- 
tion function to the ones at 70 fim, but with a temperature 
range which varies less strongly with redshift, apart from the 
very hot (A pGa k < 40 fim) SED region. On the other hand for 
the longer wavelength surveys the likelihood of detection is 
significantly more dependent on SED type. BLAST is highly 
sensitive to colder SEDs, but sources with A pea k= 60-90 fim 
can still be detected at 250 fim up to z=l. Down to 2mJy 
SCUBA should be able to detect ULIRGs over the entire 
A P cak-z parameter space, primarily because of the negative 
iv-correction which offers a large advantage above z~l. 
However above z=0.5, there is a strong predisposition to- 
wards colder (A pea k>70 fim) objects. In addition, the large 
area below A pea k=70 fim covered by the blue regions in the 
BLAST and SCUBA panels suggests that warm SEDs can 



only be recovered if they have sufficient submm flux, i.e. an 
additional cold temperature component. 

How does the dust temperature distribution of the 
BLAST, SCUBA and Spitzer samples compare in the 
context of their selection characteristics? The S09 sam- 
ple's A pea k distribution shifts from longer to shorter val- 
ues with increasing redshift, following the shape of the se- 
lection function, which shows increased sensitivity towards 
warm SEDs above z=0.7. Although above z=0.8, where 
the survey becomes more sensitive to A pca k < 70 fim, most 
ULIRGs are warm, it is interesting to note that in the unbi- 
ased region (z<0.7) most ULIRGs are cold. We calculate 
(A pea k) =99±6 fim with an intrinsic dispersion a=31 fim; 
(T) =26±2K with <t=8 K. 

The D09 sample consists of sources detected at 5a in at 
least one BLAST band. Almost all are also detected at 70 
and 160 fim, although the MIPS counterparts have fainter 
flux densities than the S09 sample. The average peak wave- 
length of the D09 sample is consistent with that of S09, 
although the characteristics of the BLAST ULIRGs cannot 
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Figure 4. SED peak wavelength (grey body temperature on the right-hand y-axis) versus total infrared luminosity for the Spitzer /70/jm 
sample of Symeonidis et al. (2009), which is at 0.4<z<1.2, with (z)=0.9. The solid and dashed lines are the L-T relation and lcr limits 
from Chapin, Hughes & Aretxaga (2009). Filled symbols indicate the true luminosities. Empty symbols correspond to the luminosities 
divided by (1+z) 3 . Here, we demostrate that evolution of the L-T relation with redshift cannot be responsible for the increased presence 
of cold ULIRGs and range of temperatures we see in the distant Universe. 



be evaluated in an unbiased part of the A poa k-z space. For 
D09 we calculate (A pca k) =89±4 with an intrinsic dispersion 
(j=15 fim; (T) =29±1 K with a=5 K. Note that the range in 
Apeak is slightly different for the MIPS and BLAST samples, 
with the former extending to lower and higher A pca k. In ad- 
dition, their redshift distributions are different, with most 
S09 ULIRGs being below z=l and most D09 ULIRGs above 
z=l. The average peak wavelength of the C08 SCUBA sam- 
ple is also in agreement with D09 and S09. For C08 we calcu- 
late (Apeak) =91±4 /im with an intrinsic dispersion a=ll /im 
and (T) =29+1 K with <r=3K, however note that (i) there 
are large upper errors on the temperature measurements, 
(ii) we are not probing an unbiased part in the selection 
function and (iii) there are additional selection effects due 
to the radio selection (see also Blain et al. 2004b; Coppin et 
al. 2009). We point out that the A poa k distribution of these 
high redshift samples extends to much longer peak wave- 
lengths than what is seen for the local ULIRGs in Fig. [21 
resulting in a higher average A pea k by at least 2 a (here a is 
the dispersion in the local sample, equal to 12 /im). 



3.2 An L-T relation at high redshift? 

In S09 we showed that a luminosity-temperature correlation 
is not apparent in the 70/zm sample, although it is clearly 
seen in the IRAS BGS. However, it has been proposed that 
the abundance of cold galaxies in the distant Universe is a 
consequence of the local L-T relation evolving with redshift 
(e.g. Chapman et al. 2002; Lewis, Chapman and Helou 2005; 
CHA09) in the same way that the knee of the infrared lumi- 
nosity function evolves by (1+z) 9 , where q~3 (e.g. Le Floc'h 
et al. 2005, Huynh et al. 2007b). This suggests that although 
more luminous galaxies would still have higher temperatures 



than their less luminous counterparts, they would be much 
colder than the objects we see locally, i.e. the ULIRG tem- 
perature distribution would shift towards longer A pca k (lower 
T) with increasing redshift. 

In Fig. |4j we examine the applicability of this hy- 
pothesis by de-evolving the luminosities of all ULIRGs 
(12<log(Li R /L ) <13) in the S09 sample by (1+z) 3 . 
CHA09 find that once the luminosities of SCUBA sources 
are de-evolved they become consistent with the local L-T 
relation. The same is not true for our lower redshift S09 
ULIRGs (0.4<z<1.2 and (z)=0.9), as most remain well out- 
side the relation, which fails to accomodate the large scatter 
in their temperatures. This result indicates that evolution 
of the local L-T relation in the ULIRG regime cannot be 
responsible for the increased presence of cold ULIRGs and 
the range of SEDs we see in the distant Universe. 



4 SELECTION OF ULIRGS IN HERSCHEL 
AND SCUBA 2 SURVEYS 

In a similar fashion to sections [2land [3~Tl we examine the po- 
tential of current surveys to measure the variation in ULIRG 
SED types by predicting the A pea k-z parameter space that 
Herschel and SCUBA 2 will cover in different bands for 
sources of log (Lm/L Q ) = 12, 12.3, 12.6 and 12.9 (Fig.[SJ). For 
the PACS sensitivity limits we take the values from Poglitsch 
et al. (2010): ~4 and lOmJy for the 70 and 160^m bands re- 
spectively for a 10' x 15'scan map with 5<r/30h integration. 
We do not include the 100 /im band, as the shape of its selec- 
tion function is similar to that at 70 and 160 /jm, effectively 
a very slow transition between the two bands. The SPIRE 
(Griffin et al. 2010) flux density limits of 20 fim for all bands, 
roughly correspond to the 3 a confusion limit (Nguyen et al. 
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Figure 5. The A pea i^-redshift parameter space for surveys with Herschel and SCUBA for log (Ltr/Lq) = 12, 12.3, 12.6 and 12.9 ULIRGs. 
For more details, see caption in figures [T] and [3] 



2010). Finally we assume a SCUBA 2 flux density limit of 
2mJy as for SCUBA in Fig. [3] 

Compared to MIPS 70 /im ultradeep data in the Great 
Observatories Origins Deep Survey (GOODS, e.g. Frayer et 
al. 2006) and other Far-Infrared Deep Extragalactic Legacy 
(FIDEL, PI: Mark Dickinson, NOAO) survey fields, obser- 
vations with the PACS 70 /jm band will mainly benefit from 
higher spatial resolution, but not from a significant increase 
in sensitivity; in terms of the latter, the superiority of PACS 
over MIPS lies largely at 100 and 160 /im. Fig. [5] shows 
that up to z=l all ULIRGs will be detected in at least 



one PACS band. On the other hand, up to the same red- 
shift, SPIRE heavily favours colder SEDs, and ULIRGs with 
A P cak>90 / iim are expected to have counterparts above the 3 a 
confusion limit in all bands. With respect to warmer SEDs, 
the extent of the blue shaded area suggests that they will 
only be detected if they have enough flux at submm wave- 
lengths, a trend which persists even for the more luminous 
(logLiR > 12.5) sources, especially at the longer wavelength 
bands. 

Between l<z<2, the PACS 70,'im selection becomes 
heavily dependent on SED shape, with many of the 
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log Lir < 12.5 objects not having 70 /im counterparts and 
a stronger tendency towards warmer SEDs for the more lu- 
minous (log Lir > 12.5) ULIRGs. The shape of the 160 /im 
selection function behaves in a similar fashion at low lu- 
minosities, however for the more luminous ULIRGs it be- 
comes relatively unbiased with respect to temperature. In 
the same redshift interval, SPIRE misses many of the lower 
luminosity sources and displays a strong tendency towards 
cold SEDs for the more luminous ULIRGs. As a result, the 
log Lir >12.5 sources should be detected in more than one 
SPIRE band if cold and in more than one PACS band if 
warm, although with certain SED shapes it is also likely 
that they will have both PACS and SPIRE counterparts. 

The picture is dramatically different at 2<z<3, where 
only the most luminous ULIRGs (logL~12.9, borderline 
HyLIRGs) will be identified in sufficient numbers. However 
their detection rate is again dependent on SED shape. For 
certain SED shapes, e.g. for sources which peak below 90 fim 
but have strong submm flux and a shallow mid-IR contin- 
uum, we can expect detections in all Herschel bands. Nev- 
ertheless, for a large fraction of ULIRGs in that redshift in- 
terval there will not be uniform and systematic information 
available on their SEDs. 

SCUBA 2 is expected to perform better than its prede- 
cessor and due to its much higher survey speed and wider 
field of view it should be able to fill out the A pC ak-redshift 
parameter space considerably better. Note that its selection 
function favours the cold ULIRG population uniformly at 
all redshifts, however its sensitivity towards warm ULIRGs 
is minimal, with a strong dependence on SED shape. The 
warm A pea k-z parameter space is covered better at high Lir, 
as the more luminous sources will have enough submm flux 
even if warm. 

It is interesting to note that at high Lir and around 
z=2, the SPIRE 250 /im selection turns over to display in- 
creased sensitivity to warmer galaxies, taking on a shape 
similar to that of the 160 /im selection function. At z ^ 2, 
250 /im corresponds to rest-frame ^ 83 fim, which means 
that the selection will henceforth favour warmer galaxies. 
This effect is only apparent in the more luminous sources, 
as lower luminosity ULIRGs cannot be detected at those 
redshifts above the 3 a confusion. A similar turnover also oc- 
curs at 160 /im, however the transition is faster in this case, 
as shorter peak wavelengths (and warmer temperatures) are 
already attainable by z=l. Note that the turnover that vari- 
ous selection functions undergo at high redshifts (and hence 
high ULIRG luminosities in a flux limited sample) can fab- 
ricate a correlation between luminosity and temperature. 
Nonetheless, this correlation is not an inherent property of 
high redshift IR populations, so caution should be taken not 
to confuse it with the intrinsic L-T relation we see in the lo- 
cal Universe (see earlier sections) — e.g. see Ivison et al. 
(2010) where sources in a BLAST 250 /im-selected sample 
need to be both more luminous and hotter in order to be 
detected at higher redshift. 



5 DISCUSSION 

We have investigated how selection effects influence our view 
of the ULIRG population both locally and at high redshift. 
Using a library of infrared SED models, we examined the se- 



lection characteristics of various infared and submm surveys 
with IRAS, Spitzer, BLAST, Herschel and SCUBA and were 
able to predict the range of ULIRG SED types detectable at 
a particular redshift, given a survey's flux density limit and 
wavelength of selection. We subsequently evaluated observa- 
tions of local and distant ULIRGs against these predictions. 
SEDs were characterised in terms of their rest-frame peak 
wavelength (equivalent to dust temperature) and total in- 
frared luminosity (in the 8-1000 /im range). A survey was 
then classed as complete and unbiased over a redshift range 
where ULIRG SEDs with peak wavelengths in the 30-150 /im 
range (T=17-87K) have equal likelihood of being detected. 

Although the sensitivity of IRAS does not enable high 
redshift observations, our analysis shows that 60 /im sur- 
veys can be unbiased over a given redshift interval, the 
size of which increases with Lir and survey depth. More 
specifically, down to a 60 /im flux density limit of 0.36 Jy, 
we find that IRAS is sensitive to all ULIRG SED types 
and luminosities in the local (z^O.l) Universe. Surveys 
with Spitzer /MIPS and Herschel /PACS which target dis- 
tant galaxies, have the advantage of remaining unbiased 
over a large fraction of the redshift range that they probe; 
on average they are complete up to z~l and their selection 
function displays little dependency on the flux density limit. 
On the other hand longer wavelength surveys with BLAST, 
Herschel/ SPIRE and SCUBA which work in the sub-mm 
regime, are significantly more sensitive to colder ULIRGs 
over the majority of their accessible redshift range. This is 
partly due to confusion in the submm limiting their depth 
and partly due to the fact that submm wavelengths probe 
emission from very cold dust at the steep Rayleigh- Jeans 
part of the SED. 

Down to the 3 a confusion limit, the synergy between 
PACS and SPIRE should enable complete characterisation 
of ULIRG SEDs up to z~l. Above z~l, the two instru- 
ments will likely sample different parts of the ULIRG pop- 
ulation but with significant overlap. For the more lumi- 
nous (logLiR >12.5) sources, surveys with Herschel and 
SCUBA 2 collectively have the potential to cover the whole 
distribution in A poak up to z~2. At redshifts beyond 2 how- 
ever, characterisation of the ULIRG population becomes 
highly incomplete, as many warm ULIRGs will not be de- 
tected and many cold ULIRGs will potentially have submm 
counterparts in fewer than 2 bands so their SEDs will not 
be measurable. 

In the context of the various surveys' selection func- 
tions we made comparisons between local and high red- 
shift ULIRGs. Using a near complete sample from the deep 
0.36 Jy 60 /im IRAS survey, we firstly examined the local 
(z^O.l) ULIRG temperature distribution. Our results show 
that, despite a uniform sensitivity to all ULIRG temper- 
atures, the average peak wavelength of local ULIRGs is 
short — (Apeak) =62 fj,m; (T) =42 K. In the same luminosity 
range (12s? log (L m /L©)<12.5), the MIPS/70 /im ULIRGs 
from Symeonidis et al. (2009) are found to be cold up to 
the redshift where the survey is unbiased (z=0.7), while 
the whole population has an average peak wavelength of 
(Apeak) =99 /im ((T) =26 K), 3<r higher than the local sam- 
ple. Although, we are not able to study BLAST (Dye et al. 
2009) and SCUBA (Coppin et al. 2008) samples in unbiased 
redshift ranges, as these are too small for current flux den- 
sity limits, we find them to have consistent temperatures to 



Selection of ULIRGs in Infrared and Submm Surveys 9 



the MIPS/70 sample — (A pcak ) =89 and 91 /an ((T) =29 
and 29 K) respectively. We note that, besides a higher value 
in the average SED peak wavelength, both the scatter and 
the range in A pca k is larger at high redshift. 

Given the differences in the A pca k distribution between 
local and high-redshift ULIRGs, one might also ask whether 
the correlation between IR luminosity and dust tempera- 
ture, which has been observed in the local Universe, is also 
applicable at high redshift. The local L-T relation dictates 
that in the ULIRG regime the average temperature should 
be high: the L-T relation in Chapin, Hughes & Aretxaga 
(2009) extrapolated to ULIRGs predicts temperatures in the 
40-65 K range within 1 a. Although our results regarding the 
local 0.36 Jy sample are in agreement with this temperature 
range, the extrapolation in L-T for the ULIRG regime fails 
to accomodate both the large scatter in temperature and 
the low average temperature seen at high redshift: (i) only 
16 per cent of the 70 /im sample lie within the local L-T 
relation 1 a limits, whereas no sources in the SCUBA and 
BLAST samples lie within the local L-T relation and (ii) the 
3 samples display no apparent connection between Lir and 
Apeak- Including results presented here and in S09, ULIRGs 
at z ~ 1 and beyond which span a large range in dust temper- 
ature have been discovered in many surveys. Cold ULIRGs 
have been identified with ISO (e.g. Chapman et al. 2002), 
Spitzer (e.g. Kartaltepe et al. 2010) and SCUBA (e.g. Chap- 
man et al. 2005). Warm ULIRGs, many of which do not have 
appreciable submm detections, but are detected at shorter 
infrared wavelengths or other parts of the spectrum alto- 
gether (such as radio or optical), have also been shown to 
exist (e.g. Blain et al. 2004a; Chapman et al. 2004; Buss- 
mann et al. 2009, Casey et al. 2009, Younger et al. 2009). 
The collective results from these studies establish a large 
range in ULIRG temperatures, suggesting that they do not 
conform to a universal luminosity-temperature correlation. 
Note that although such a correlation is sometimes observed 
(e.g. Blain, Barnard & Chapman 2003; Chapman et al. 2004, 
2005; Ivison et al. 2010), it does not have a physical origin. 
Rather, the surveys' selection functions imply that at above 
a certain redshift, ULIRGs will need to be both more lu- 
minous and hotter in order to be detected. In theory, the 
selection function of all surveys should turn to favour warm 
sources at some redshift, however the longer the wavelength 
of selection, the higher is the redshift that this will occur at. 

In order to explain the abundance of cold ULIRGs at 
high redshift and the lack of them in the local Universe, it 
has been suggested that the local L-T relation evolves in the 
same way as the knee of the luminosity function, by (l+z)' 5 , 
i.e. galaxies become colder at high redshift. However, the 
large observed scatter in ULIRG temperatures and the lack 
of an observable trend between luminosity and temperature 
amongst distant ULIRGs is a strong argument against this 
proposal; we demonstrate this in Fig. [4] and Section 13.21 
So how can one characterise the discrepancy in the tem- 
perature distributions of local and high redshift samples? A 
recent study by Seymour et al. (2010) examined the evo- 
lution in infrared luminosity density (IRLD) with redshift 
as a function of dust temperature, by splitting sources into 
warm and cold (at A pca k=90 (im) and examining their contri- 
bution to the IRLD separately. Seymour et al. showed that 
cold galaxies undergo stronger evolution and are primarily 
responsible for the increase in infrared luminosity density 



with redshift, a result which can describe the differences in 
the ULIRG temperature distributions we see locally and at 
high redshift. If cold ULIRGs evolve faster than their warm 
counterparts, then we would expect their relative numbers 
to diminish with decreasing redshift. Although, the division 
of galaxies into two groups, warm and cold, is possibly too 
simplistic, this picture does highlight a possible direction for 
future work on this topic. With clean samples from carefully 
tailored surveys over redshift ranges where we have shown 
them to be unbiased, Herschel should be key in fully charac- 
terising ULIRG SEDs and hence investigating the evolution 
of the IR luminosity function as a function of temperature. 



6 CONCLUSIONS 

Our conclusions are summarised below: 

(i) All surveys, independent of the wavelength of selec- 
tion, have the potential of being complete and unbiased with 
respect to ULIRG SEDs over a certain redshift range, the ex- 
tent of which increases with survey depth and ULIRG lumi- 
nosity. As it stands, short wavelength (A <200 (im) surveys 
with IRAS, Spitzer /MIPS and Herschel /PACS are sensitive 
to all SED types in a large temperature interval (17-87 K) 
over a substantial fraction of their accessible redshift range. 
On the other hand, long wavelength (A >200 /im) surveys 
with BLAST, tfer.se/id/SPIRE and SCUBA are significantly 
more sensitive to cold ULIRGs, disfavouring warmer SEDs 
even at low redshifts, partly due to confusion in the submm 
limiting their depth and partly due to the fact that submm 
wavelengths probe emission from very cold dust at the steep 
Rayleigh- Jeans part of the SED. 

(ii) By examining a near complete IRAS 60 /im sample 
up to z=0.1, where the survey is sensitive to all SED types 
with average temperatures in the 17-87 K range, we confirm 
that cold ULIRGs are rare in the local Universe. The local 
ULIRG SEDs have an average peak wavelength of 62±2 fim 
with an intrinsic dispersion <r=12/xm which translates to a 
greybody temperature of 42±2 K with a=8 K. In compari- 
son, for higher redshift Spitzer /MIPS, BLAST and SCUBA 
samples, we calculate (A pea k) =99±6 fim with an intrinsic 
dispersion a=31 /im, 86±4 with o = 15 /jm and 91±4 /im with 
a=ll ^m ((T)=26±2K with (7=8 K, 29±1K with a=5K 
and 29±1K with o=3K) respectively, indicative of an in- 
crease in the relative number of cold galaxies at high red- 
shift. 

(iii) The local L-T relation which describes luminous ob- 
jects as warmer than their less luminous counterparts, does 
not apply to ULIRGs in the distant Universe, as the scat- 
ter in the temperature distribution of high redshift ULIRGs 
is much higher than the ULIRG temperatures predicted by 
the L-T relation. Only 16 per cent of the MIPS sample falls 
within 1 cr of the L-T relation, compared to 70 per cent of 
the local ULIRG sample. 

(iv) The discrepancies between local and high redshift 
samples, cannot be reconciled by assuming that the local L- 
T relation evolves with redshift in the same way as the knee 
of the IR luminosity function, i.e. by (1+z) 3 . Instead, we find 
that these differences can adequately be described by the 
Seymour et al. (2010) results, which show that evolution of 
the infrared luminosity function is temperature dependent, 
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with cold galaxies evolving at a faster rate, implying that 
their relative number will decrease with decreasing redshift. 
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Figure Al. SED peak wavelength versus rest-frame colour C 
[=log Leo/Lioo where L A is in W/Hz] for the 27 local IIFSCz 
ULIRGs. The solid and dotted lines are equation IA1I (and corre- 
sponding errors). 



APPENDIX A: THE RELATION BETWEEN 
COLOUR AND SED PEAK WAVELENGTH 

In order to use the CHA09 L-T relation, which uses colour, 
instead of SED peak wavelength, we establish a conversion 
between the two quantities. Using our fitted SED templates, 
we calculate rest-frame Leo/Lioo colour (C) for the 27 local 
IIFSCz ULIRGs (where h\ is in W/Hz for direct comparison 
with the C [= log (feo/fioo) where f\ is in Jy] definition used 
in CHA09). As expected, C and A poa k follow a one to one 
relationship which can be fit with the following equation 
(using a biweight estimator): 

Ap C ak = 58.25(±1.35) - 72.82(±12.2) C (Al) 

where C=log (Leo/Lioo), shown in figure |A"T1 



